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abstract 

An  antisy^aetrically  laminated  angle-ply  plate  is  optimized  with  the 
objectives  Of  maximizing  the  fundmaental  eigen  frequency  and  the  distan, 
between  two  consecutive  natural  frequencies.  The  fonnulation  includes 
the  contribution  of  the  shear  deformation,  but  neglects  the  in-pl,ne  an 
rotary  inertias.  The  design  variables  are  the  fiber  orientations  and 
the  thicknesses  of  individual  layers.  The  design  problems  are  cast 
Into  a  mathematical  progranming  fomat  and  solved  by  using  a  quasi- 
Newton  function  maximization  algorithm.  A  penalty  function  method  is 
employed  to  maximize  the  fundamental  frequency,  subject  to  lower  bound 
constraints  on  higher  order  frequencies.  Numerical  results  are  presente 
for  laminates  constructed  of  high  modulus  fibr^e  reinforced  materials 
«.d  the  effects  of  various  problem  parameters  on  the  efficiency  of  the 
designs  are  Investigated,  it  is  shown  that  the  design  variables  may 
not  be  determined  optimally  i  f  the  effect  of  shear  deformation  is 
nogldcted.  Moreover,  it  wa,  observed  that  the  classical  plate  theory 
loads  to  erroneous  results  in  optimal  material  selection  problems. 
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1.  IHTRODUCTlOff 

Ir.  a  recent  article  Bert  and  Chen  til  studied  the  effect  of  shear  defor= 
mation  on  the  natural  frequencies  of  an  antisymnetric  angle-ply  laminate 
and  showed  that  the  classical  plate  theory  (CPT)  tends  to  overestimate 
the  values  of  frequencies.  Specifically,  rectangular  plates  on  simple 
supports  were  considered,  and  the  effects  of  in-plane  and  rotary  inertias 
were  included  in  the  formulation  til.  Subsequently,  the  same  plate 
problem  was  also  analyzed  by  means  of  a  finite  element  method  [21.  in 
the  present  article,  a  rnanber  of  optimal  design  problems  for  similar 
structures  are  solved  by  employing  the  shear  deformation  theory  (SDT) 
given  in  [IJ.  He  neglect  the  effects  of  In-plane  and  rotary  Inertias, 
which  were  shown  to  have  little  effect  on  the  fundamental  frequency  even 
In  the  case  of  relatively  thick  plates  [11.  In  particular,  we  consider 
the  problems  of  maximizing  the  fundamental  frequencies  of  angle-ply 
laminates  with  and  without  lower  bounds  on  higher  order  frequencies  and 
of  maximizing  the  distance  between  two  consecutive  frequencies. 

The  design  variables  are  taken  as  the  fiber  orientations  and  thicknesses 
Of  individual  layers.  Optimization  problems  are  foranjlated  as  mathematical 
progra»«.lng  problems  and  solved  by  using  a  quasi-Newton  nonlinear  function 
maximization  algorithm.  In  the  constrained  optimization  problem,  the 
lower  bounds  on  higher  order  frequencies  are  Incorporated  into  the 
formulation  by  means  of  a  penalty  function  technique. 

Optimal  designs  of  laminated  plates  were  given  in  several  studies  13-81 
with  respect  to  natural  frequencies  in  which  the  effect  of  shear  deforma^ 
tion  was  neglected.  He  refer  the  reader  to  [7]  for  works  concerning  the 
optimization  of  laminates  with  respect  to  oUier  objectives,  such 
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for  the  maximum  separation  of  frequencies  were  obtained  for  isotropic 
bea^s  by  Olhoff  tl/.IS]  by  .aximi.ing  the  higher  order  frequencies. 
the  case  of  composite  structures,  there  does  not  seen,  to  exist  any  study 
on  this  problem.  The  present  article  gives  designs  for  maxim*,  distance 
between  the  lowest  and  the  next  lowest  two  natural  frequencies.  In 
115.17,18].  it  was  observed  that  the  higher  order  frequencies  approach 
each  other  when  the  structure  is  designed  with  respect  to  frequency 
separation.  The  same  phenomenon  is  also  observed  in  the  present  case  and 
seems  to  be  a  general  characteristics  of  such  designs.  We  offer  an 

explanation  for  this  in  terms  of  the  topology  of  the  frequency  surfaces 
in  the  design  space. 

2.  PROBLEM  FORHULATIQW 

We  consider  a  simply  supported  rectangular  plate  of  length  a.  width  b. 
thickness  h  and  mass  density  p.  The  plate  is  composed  of  an  even  number 
of  orthotropic  layers  of  thickness  H^.  the  fibers  of  which  are  oriented 
alternately  at  angles  ano  -e^  and  are  placed  anti syametri cal ly  with 
respect  to  the  middle  surface.  Plates  with  these  characteristics  are 
conmonly  known  as  antisyimetric  angle-ply  laminates  tl.2.7.8].  The 
equations  of  motion,  which  include  the  effect  of  shear  defcnnation.  were 
given  in  tU  for  freely  vibrating  angle-ply  laminates.  We  adopt  the  same 
fonaulation  in  the  present  study  but  neglect  the  effects  of  in-plane  and 
rotary  inertias.  The  equations  are  presented  in  the  /Ippendix.  We  obtain 
a  non-dimensional  fom  of  these  equations  by  introducing  the  following 
dimensionless  quantities: 


X  =  x/a.  y  =  Y/b.  u  =  U/a.  v  *  V/b.  w  =  W/a 
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r  =  a/b,  p  .  a/h,  hj^  =  Hj^/h, 


(1) 


vrfiere  i  and  J  are  integers  and  k  refers  to  the  layer  number.  In  equation 


(1),  X  and  Y  are  Cartesian  coordinates  parallel  to  the  respective  plate 
edges;  U,  V,  U  are  the  displacement  components  in  the  directions  of  the 
XYZ  system;  A.j,  and  O.j  are  the  laminate  stiffnesses  given  by 

h/2 

<*ij*®ij*'’ij>  =  «.  (2) 

where  denotes  the  plane  stress  reduced  stiffness  components  of  the 
k-th  layer  given  In  the  Appendix.  is  the  transverse  modulus  of 
elasticity,  but  could  be  replaced  by  any  reference  modulus.  Introducing 
(1)  into  (Al),  we  obtain  the  following  non-dimensional  system  of  coupled 
partial  differential  equations  governing  the  free  vibration  of  the  shear 
deformable  laminate: 


+  »66'""''yy  "  ^  •’iBP'V'fxx 

"  ‘•zr/VNyy  "  ♦xy  =  0 

<*12**66>“xy  ♦  »6f,'-'^''xx  *  *22^  ♦  2‘>26P'‘  ^xy 

*  '’l6''''p'\x  ♦  Vp’’  V  '  0 

■'s  ^ss^xx  ♦  •‘u  *  K  a.,,/  fy 


(3) 


(4) 


(5) 
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*  2^eP'‘»xy  - 

*  ^66P  ’’xx  *  *'22'’  P  ^  ’^yy  *  '‘ii®M<i  ’’  (®) 

*  <''l2*'*66>'1^’^  **y  =  » 

2^6’T>‘'“xy  ♦  '',6'''’P*‘''XX  ^  '>26'T>‘‘v 

■'‘^55*'x  *  {‘*12*‘’66>'T’'^  ^y  ♦  <*iiP*^  ♦xx  (2) 

*  "ee-^'p'^  V  -  ♦  =  0 

where  ♦  and  V  are  the  slopes  in  xz  anJ  yz  planes  respectively;  and  k^ 
are  shear  correction  coefficients  and  t  denotes  tins. 

Me  consider  the  same  simply  supported  plate  boundary  conditions  as  the 
ones  given  in  tH.  In  terms  of  dimensionless  quantities,  these  are 

u(0.y)  vi(l.y)  =  0.  Hp(x,0)  =  Ng(x.l)  =0 

NfifOry)  '  “O’  ’'i"'®)  =  vCx.l)  =  0 

w(0,y)  =  w(l,y)  =  w(x,0)  =  w{x,i)  =  0 
Hi(0.y)  =  Hj(l,y)  =  M,^(x.O)  =  M^(x,l)  .  0 
f(O.y)  '  T(l,y)  =  0,  ♦(x.O)  =  ♦(x.l)  =  0 

where  the  dimensionless  stress  couples  Mj  and  and  in-plane  stress 
resultant  are  given  by 

"c  '  *66<'’*‘''x^™y)  ^  bleP'*  *X  *  '>26''P'’  ’’y 

”1  =  •'i6('’'’vv  r '’1/'  %  ♦  '*i2'-p’'  \  (3) 

”2  -  V'’VV  "  '*12P'‘  ^  »  ''22'T’‘’  -fy 

The  objectives  of  the  study  are  to  maximize  the  fundamental  eigenfreqiiency 
Oj  and  the  distance  between  the  consecutive  natural  frequencies  of  the 


.rrr  ’’*  ”'  "■  *■''’  '”"•  “■  ""••  •“  f- 
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»  •».,..  for  ™„™, 
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problem 

*MX  Q  .  |(=1  ?  i/yo 

ti  * ^.^-...K/e, 

®k  ‘^k 

where  the  design  variables  0^  and  h.  satisfy 

K 

>0.k=lP  if 

k.l  ^  k  1,2*.. ..K,  (IPJ 

with  K  denoting  the  total  nuiaher  of  layers. 

we  note  that  doe  to  antisy^^try.  we  have  e,  = 

s  a  result  of  these  require, nents  and  (12).  „e  have,  in  effect  t:/2 

vanabies  of  n,  and  (K/2,-1  variables  of  h,.  Thus,  the  total  „ua.er  of 
effective  design  variables  is  K-i. 
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ConstralnelDe^^  Oeten.,'™,  the  solution  of  the  -«xi™it,tlon 
problem 


subject  to 


.  k=l.2.....K/2 
”l,  'i  »  2.  . 


(13) 


(M) 


where  end  h,  satisfy  (11)  and  (12)  respectively  and  denotes  the 

1,-th  order  frequency,  H  a  specified  lower  bound. 

Mfgn  for  waxlmum  frequenry^eea!^;  Oetennine  the  solution  of 
tnsxfnif zdtion  problem 


®k  ’*’k 

where  0^  and  satisfy  (II)  and  (12).  respectively,  and  n,  denotes  the 
1-th  order  frequency. 

..  *»n, 

6j^  and  bj^. 

3-  HETIIOO  of  SOLUTIOH 

The  foMowin,  set  of  displacement  and  rotation  functions  satisfies  the 
differential  equations  (3)-(7)  and  the  bo.«.dary  condition,  (8)  ,  ,1: 

u  =  Sin  Hhfx  COS  nny  COS  n  t 

mn 

^  cos  rmix  sin  nny  cos  I2  t 

mn 

^  cos  n  t 

nin 

*  ’’’mn  cos  ntiy  cos  ft  t 

nin 

^  ~  ^ran  cos  ft  t 

mn 


> 

f 


a 


where  the  frequency  of  t>»e 


eigentnode  (in,n). 


'"^ertlon  of  (16,  i„to  (3,.(;,  ,,,,3 
•Jgebrale  eiuatfons  of  tho  form; 


set.  of  linear  homogeneous 


where  S  =*  t  u  y  w  ♦  *  -.t  , 
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M  •»!< 

Sb  '  * 

Sr,  '  - -nP-V  -  d^^rVV  - 

f!7  V  -''-^‘onless  elgenfro^ooncy 

‘  mn  •’>  ‘be  formula  '  mn  '^'ates 


"■mn  *PaV/¥' 

A  nontrivial  solution  of  (17,  exfsu  If 

<tet  C  =  0, 

Which  leads  to  explicit  expression  for  „  . 

mn  *  * ' * 


(19, 


(20, 


.» 

f 
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“4  <'‘^55“(''35^''33>  ♦  •'^M/SCFaM+rPFjj))  (21) 

where  F,j  denotes  the  cofactor  of  the  eleoent  C,,.  Ue  note  that  an  explicit 
expression  for  can  be  obtained  due  to  the  fact  that  the  In-plane 
and  rotary  Inertias  are  neglected.  Otherwise  the  determlnental  equa'tlon 
(20)  becomes  a  fifth  order  polynomial  in  -uni- 


Maximization  of  the  expression  (21)  with  »=n*l  and  with  respect  to  the 
design  variables  0^  and  which  are  subject  to  the  constraints  (11) 
and  (12),  leads  to  the  solution  of  ttw  unconstrained  design  problem  (10). 
In  the  case  of  the  constrained  design,  the  fundamental  frequency  is 
to  be  maximized  subject  to  constraints  w.  2r  Wp  ,  t.  >  2,  1«1,2,,..  on 
f|-th  order  frequencies  This  type  of  problem  was  already  treated 

by  Adali  f8,191  in  the  context  of  laminated  plates  and  vibrating  beams 
by  employing  a  penalty  function  method.  The  procedure  converts  the  con^ 
strained  optimization  problem  into  an  unconstrained  one  by  treating  the 
prograBHiing  problem 


max 

0..  h.  ” 


■  1 

+  J:  —  max  fO,  6). 

1=1  ^1  S 


1)  ♦ 
i 


(22) 


where  >  0  are  specified  small  constants.  For  further  details  on  the 
penalty  function  method  we  refer  the  reader  to  (19).  Solution  of  problem 
(22),  subject  to  (II)  and  (12)  provides  the  constrained  optima?  design 
of  the  laminate. 


We  note  that  both  the  constrained  design  problem  and  the  maximum  frequency 
separation  problem  involve  the  evaluation  of  higher  order  frequencies. 
Although  the  fundamental  frequency  is  given  by  for  all  the  problem 
parameters,  the  modal  wave  numbers  ra  and  n  depend  on  the  specific  problem 
parameters  in  the  rase  of  higher  order  frequencies.  We  compute  the  t-th 
order  frequency  from 

o-j  .  Bln  subject  to  ,  (23, 
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which  constitutes  a  discrete  minia,1„tion  probiern  over  the  integers  n,  and 
n. 

^  DESIGN  FOR  MAXIMUM  FUNDWENTAL  FREQUCHry 

In  this  section,  numerical  results  for  the  optical  design  of  an  angle-ply 
I».nate  for  maximum  are  given,  and  the  effect  of  various  problem 
parameters  on  the  efficiency  of  a  design  is  investigated.  «e  treat  two 
different  kinds  of  graphite  epovy  plastic  with  material  constants  given 


Material  I:  E^/e^  .  40.O.  G^^/E^  =0.6 

Gj^/Et  =  0.5.  .  0.25 

and  '** 

Material  II:  E^/E^  .  iS.O.  G^^/Ej  =  0.5 
Gt^/Fj  =  0.2.  .  0.25. 

The  shear  coefficients  are  taken  as  k,^  =  k^  ,  S/6. 

•lency  of  an  optimal  design  is  assessed  by  comparing  it  with  the 
corresponding  standard  plate,  which  we  define  to  be  composed  of  layers 
Of  egual  thickness  with  the  fibers  oriented  alternately  at  45  and  -45 

degrees.  The  efficiency  index  is  defined  as 

fff  '  100({.,.^p/,„j)-  I). 

Which  gives  the  percent  increase  in  the  fundan«ntal  freguency  of  the 

optimal  plate  as  compared  to  the  fundamental  frequency  of  th^standard 
plate.  ^ 

r.,.l  ... 

f.r  „  . . .  ^ 


those  for  a/h  =  «  are  obtained  by  employing  the  classical  plate  theory. 

We  observe  that  the  effect  of  shear  deformation  Is  to  decrease  the  effi¬ 
ciency  of  the  design  and  this  effect  becomes  more  pronounced  as  the  a/h 
ratio  becomes  smaller.  Exactly  the  same  situation  was  observed  in  the 
optimal  thickness  design  of  one-dimensional  structures  that  included  shear 
deformation  [9-121. 

Fig. 2  gives  the  efficiency  curves  for  laminates  made  of  material  II  with 
the  same  problem  parameters  as  in  Fig.l.  A  comparison  of  Figs. I  and  2 
indicates  that  material  I  produces  less  efficient  designs  for  a/h  <  10 
when  a/b  <  I  and  for  a/h  <  40  when  a/b  >  1  than  material  II,  but  this  is 
not  the  case  for  a/h  =  «.  An  important  implication  of  this  observation 
is  that  an  assessment  of  efficiency  of  plates  constructed  from  different 
materials  would  be  completely  misleading  if  the  shear  deformation  were 
neglected. 

Tables  1  and  2  provide  the  values  of  fundamental  frequencies  and  fiber 
orientations  of  optimally  designed  laminates  made  of  materials  I  and  II 
respectively.  We  note  that  the  results  for  optimum  0. 's  are  given  for 
only  one-half  of  the  laminates  because  of  tlte  anti  symmetry.  Table  1  in* 
di cates  tiiat  the  fiber  orientations  of  optimal  plates  depend  on  the  a/h 
ratio,  and  consequently  SDT  and  CPT  yield  different  results.  We  observe 
from  T«:ble  2  that  these  differences  are  less  pronounced  for  material  II. 
Thus  the  extent  of  the  difference  in  optimum  0^*s  obtained  by  SOP  and  CPT 
is  closely  related  to  the  material  properties  and  mostly  to  E|^/Fj  ratio. 

Next,  we  investigate  the  amount  of  decrease  in  the  maxintum  fundamental 
frequencies  of  optimal  designs  as  a  result  of  taking  shear  deformation 
into  account.  Let  and  denote  the  fundamental  frecu'ncies  of  op® 
timal  laminates  obtained  by  SDT  and  CPT  respectively.  Then,  the  percent 
decrease  in  efficiency  as  a  result  of  employing  SOI  rather  than  CPT  is 


given  by 


“eff  '  ■  '"SD  '  "Cp)-  (25) 

The  curves  of  0^^^  plotted  against  the  aspect  ratio  are  given  in  Fig.3 
for  four-layered  laminates  optimized  with  respect  to  e..  We  observe  that 
the  efficiencies  obtained  from  CPT  become  less  accurate  as  a/h  decreases 
and/or  a/b  increases.  Moreover,  the  material  properties  have  relatively 
small  effect  on  the  final  result. 

In  the  above  results,  the  laminates  were  optimized  with  respect  to  fiber 
orientations  only.  The  effect  of  including  the  layer  thicknesses  among 
the  design  variable,  is  studied  in  Fig.4.  which  shows  the  differences  in 
the  efficiencies  of  optimal  plates  with  design  variables  {e,.h,)  and  (o.). 
In  the  particular  case  of  four- layered  plates  considered  In  Fig. 4.  the  ' 
specific  design  variables  are  (o^.e^.h^)  and  (O^.o^).  The  percent  differ 
rences  in  the  efficiencies  are  given  by 

““‘■''op’  - '4p’)/-s. 

where  a„d  .g)  denote  fundamental  frequencies  of  optimal  laminates 

with  three  and  two  design  variables  respectively.  We  observe  that  the 
optimum  layer  thickness  becomes  most  effective  when  aspect  ratios  an* 
close  to  unity,  and  moreover  the  increase  in  the  efficiency  depends  heavily 
on  a/h  ratio.  Another  interesting  conclusion  is  that  the  optimum  deter= 
mination  of  the  layer  thickness  increases  the  efficiency  of  a  laminate  for 

aspect  ratios  where  CPT  would  indicate  no  change.  More  specifically  CPT 

sht«s  no  increase  in  the  efficiency  of  a  design  as  a  result  of  including 
the  layer  thickness  among  the  design  variables  when  0.55  >  a/b  >  1.7. 

Table  3  gives  the  values  of  the  maximum  and  optimum  9, .0j,h^(=z{3)=-z(I)) 
for  various  a/h  and  a/b  ratios  and  for  laminates  made  of  material  I  A 
Ci^npdrison  bf  tables  i  and  3  indicates  that  optimum  values  of  o,  depend  on 
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the  number  of  design  variables. 

Fig.5  gives  the  curves  of  efficiency  plotted  against  a/h  for  laminates 
with  two.  four  and  six  layers.  He  observe  that  the  efficiency  drops  as 
the  number  of  layers  increases.  Horeover.  the  efficiencies  of  Materials 
I  and  II  «y  be  greater  or  less  than  one  another  depending  on  a/h  as 
illustrated  in  the  case  of  four-layered  laminates.  Thus  we  again  come 
to  the  conclusion  that  a  comparison  of  the  efficiencies  of  different  ma. 

terials  by  CPT  alone  is  likely  to  give  incorrect  results  for  thick  lami= 
nates. 

The  effect  of  the  modulus  ratio  on  the  efficiency  is  specifically 
Investigated  in  Fig.6.  where  the  efficiency  curves  are  plotted  against 
E^/E^  for  laminates  optimized  with  respect  to  9,  and  with  .  0.6, 

Gt^/Et  =  0.5.  ;  0.25,  r=2  and  number  of  layers  4.  He  observe  that 

the  efficiency  decreases  after  a  certain  value  of  E^/E^  as  E^/E^  increases 
when  a/h  <  40.  On  the  other  hand  CPT  indicates  a  steady  increase  in 
efficiency  for  the  range  of  values  of  E^/E^  given  in  Fig.6.  Thus,  even  a 
qualitative  estimate  based  on  CPT  as  to  the  efficiency  of  laminates  made 
of  different  materials  will  be  incorrect.  From  Fig.6.  it  appears  that 
the  effect  of  shear  deformation  increases  as  F^/E/  inc leases,  causing  the 
efficiency  to  drop  after  a  certain  point.  The  initial  increase  in  effi. 
ciency  apparently  occurs  when  the  deformation  due  to  bending  only  dominates 
the  contribution  of  the  shear  deformation. 

COWSTRAiHED  OPTlHimiOil  AMD  OPTIMAL  FREQUENCY  SEPflRATlOM 

In  many  applications,  the  maximization  of  the  fundamental  frequency  may 
be  subjected  to  additional  design  requirements  in  the  form  of  constraints 
on  higher  order  frequencies.  Usually  one  or  more  of  the  higher  order  fre= 
quencies  are  required  to  be  greater  than  certain  values.  Another  class 
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of  problems  is  the  design  of  a  laminate  that  has  maxiimim  separation  between 
the  consecutive  frequencies.  The  optimization  technique  used  in  the 
present  study  can  be  easily  extended  to  solve  both  problems. 

Table  4  gives  the  values  of  the  maximum  fundamental  frequencies  with  second 
and  third  order  frequencies  subjected  to  lower  bounds  and  respec= 
tively.  The  results  are  given  for  four- layered  laminates  optimized  with 
respect  to  only,  with  a/b  =  10.0  and  a/b  =  1,2.  When  =  24.0  and 
W.J  *  with  a/b  ‘  I,  the  laminate  has  a  double  eigenvalue  as  its  second 
order  frequency.  This  situation  changes  when  we  set  =  37.0,  and  niore= 
over  the  mode  of  the  third  order  frequency  changes  from  (2,1)  to  (1,3) 
as  w  is  Increased  to  =  38.0.  The  results  for  r=2  indicate  that  w 

J  3  2 

and  w  can  be  increased  considerably  with  little  effect  on  w  , 

11 

Numerical  results  for  the  optimal  frequency  separation  problem  are  given 
in  Table  5  for  similar  parameters.  We  observe  that  the  maximizations  of 
with  r=l,2  lead  to  double  eigenvalues  for  the 
second  and  third  order  frequencies  respectively.  This  may  be  attritn»ted 
to  the  fact  that  a  frequency  surface  in  the  design  space  usually  has  a 
local  maximum  at  those  points  where  the  surfaces  of  different  eigenmodes 
cross  each  other.  Obviously  these  points  are  also  those  where  double 
eigenvalues  occur. 

6.  CONCIUSIOMS 

The  main  effect  of  the  shear  deformation  on  an  optimally  designed  laminate 
is  to  reduce  its  efficiency.  We  note  that  the  efficiency  is  defined  in 
comparison  to  the  corresponding  antisymmetric  angle-ply  laminate  that  has 
equal  thickness  layers  with  fibers  oriented  at  alternating  angles  of  45 
degrees.  The  efficiency  decreases  with  decreasing  side- to-thickness  ratio, 
as  is  the  case  in  the  optimal  thickness  design  of  isotropic  structures  r9-12l. 
The  aspect  ratio,  number  of  layers  and  the  material  properties  have  considerabl 
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effect  on  the  efficfency  of  a  design,  as  depicted  in  Figs.  I.2.5  and  6. 

comparison  of  results  obtained  by  SOT  and  CPT  yields  interesting  insights 
Into  the  ways  the  shear  defonnation  affects  the  optimal  laminates  made 
of  advanced  filamentary  composite  materials.  For  example,  SOT  and  CPT 
produce  different  values  for  optimum  fiber  orientations  and  layer  thick, 
nesses.  Thus,  a  design  based  on  CPT  alone  could,  in  fact,  be  only  subop. 
tioal.  Deviation  from  the  optimum  increases  with  decreasing  side- to- thick 
ness  and  increasing  E^/E^  ratios  when  CPT  is  expioyed. 

Furthermore.  CPT  gives  qualitatively  Incorrect  results  when  a  cmnparison 
of  different  composite  materials  is  made.  Indeed.  CPT  Indicates  that  the 
efficiency  would  increase  with  increasing  E^/E,  ratio.  Therefore,  the 
material  with  the  highest  E^/E,  yields  the  most  efficient  design  accnr= 
ding  to  CPT.  In  fact,  the  efficiency  decreases  after  a  certain  point  with 
Increasing  E^/E^  when  the  shear  is  taken  into  account.  Thus.  SOT  indicates 
an  optimal  value  for  the  maximum  efficiency. 

The  efficiency  can  be  increased  by  including  the  layer  thicknesses  a»,na 
the  fiber  orientations  as  additional  design  variables.  This  increase  Is 
larger  when  the  side-to-thickness  ratio  is  high  and  when  the  aspect  ratio 
Is  close  to  unity,  but  it  becomes  negligible  if  the  aspect  ratio  is  too 
small  or  too  large.  CPT  yields  misleading  results  on  this  aspect  by  indi. 
eating  too  narrow  a  range  for  the  aspect  ratio  around  unity  where  optimizing 
with  respect  to  layer  thicknesses  can  be  effective. 

Constrained  designs  of  the  laminates  were  obtained  when  lower  bounds  were 
Imposed  on  higher  order  natural  frequencies.  It  was  observed  that  the 
imposed  bounds  effected  the  eigenmode  of  the  higher  order  frequencies 
at  the  optimum  point  (Table  4). 

Designs  for  optimal  frequency  separation  lead  to  double  eigenvalues  for 


:> 
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I  higher  order  frequencies  in  many  cases  (Table  5).  This  seems  to  be  due 

I;  . 

f  to  local  maxima,  which  occur  at  the  cross  sections  of  surfaces  formed 

i.  ■ 

V  by  different  eigenmodes  in  the  design  space, 

t 
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appendix 

The  equations  of  motion  governing  the  free  vibrations  of  an  antisy«etric 
angle-ply  laminated  plate  are  given  by  [11: 

t  0  .  0, 

where  p  =  fU  V  W  hr  is  the  transpose  of  the  displacement  vector  and 

£  '  ri^jJ  is  the  5  by  5  symaetric  differential  operator  matrix,  the  cc»n= 
ponents  of  which  are  given  by 

'■It  *  *11  “x  *  *66  ‘•12  °  (*12  + 

••13  '  O-  ■  (Bj5/h)0*  +  (Bj6/h)Dj 

•■IS  '  <^®16''‘')“x’’y'  ‘•22  '  *66®x  ^  *22  “v 
^23  =  0-  '•2,  =  (2B,6^h)Dx'’Y. 

‘•33  '  *''5  *55  ®X  ■  *^4  *i)i|  bJ  +  p  h  0^ 

'■34  '  '('‘4  ^q/hlDy.  Lj5  =  -(k*  Ajj/h)Dj 

‘•44  ^  (‘'66'''^)‘>X  ♦  (‘>22/''")‘>?  -  ‘'J  A4,./h’ 

*•45  '  (vOj,  ♦ 

••SS  ’  <‘’ll'''^)‘’x  *  <®66^'»^)B?  - 
where  Ojj  ?  a/ax,  =  a/aY,  0^  =  a/at. 

Plane  stress  reduced  stiffness  components  of  the  k-th  layer  are  given  by 
0{J>  =  0,1  c"  a  2q^,_.„2c2  *  j4 

•Jia’  =  (•’,,♦  0„  -  4Q^^,)sV  t  Qj^(s“  .  c'-) 

'’22’*‘»11*-^2(Q,2^2Q66)=V-^  022  ^‘‘  (A2) 

(k) 

V  •  «ir  0,2  -  V  -  Q„  t  2q^Js3c 

(k) 

V  “  {9,1  -  Q22  -  2066^'^  t  (0,,  -  4  2Q^^)S 
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c"  +  o„  s" 

•  <''it*'»22-2'>i2-2Qfie)s"c^0,g(sV). 
where  c  =  cos  0^  .  s  =  s1„ 

•’ll  '  Ojj'  '’ltV(>-'’lt'’tlJ 

<>2?  =  Q,,^  = 

%S  *  "^68  °  ®LT»  '  ''lt  ^T^^l- 

2^^  h  O-te  Voo„,.s  in  tHe  ,o„.n„d1na,  and  traosverse 

.  '“■  'u  ’> »« 

Strain  under  longitudinal  strp«;c«  r  j  « 

nessch  "U  *"«<  G^z  denote  in-p,ane  and  thici 

ness  shear  moduli  respectively. 
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TABLE  1:  Mixlnw  fundaBental  eigpnfrequencles  and  optimum  fiber  orienta- 
lions  in  degrees  of  a  four-layered  angle-ply  plate  made  of 
material  I.  with  .  0.25.  k.1.2,3,4. 


a 

H 

5 

10 

20 

40 

a 

F 

0,0 

9.46 

(0.0)* 

14.00 

(0.0) 

16.71 

(0.0) 

17.67 

(0.0) 

18.03 

(0.0) 

0,5 

9.73 

(14.4/ 

-30.0) 

14.16 

(5.2/ 

-15.8) 

16.83 

(0.0) 

17.85 

(0.0) 

18.21 

(0.0) 

1.0 

12.67 

(45.0/ 

-45.0) 

18.53 

(45.0/ 

-45.0) 

21.91 

(45.0/ 

-45.0) 

23.09 

(45.0/ 

-45.0) 

23.53 

(45.0/ 

-45.0) 

2,0 

22.55 

(68.6/ 

-61.4) 

38.93 

(75.6/ 

-60.0) 

56.62 

(84.8/ 

-74.2) 

67.50 

(90.0) 

72.84 

(90,0) 

4.0 

44.57 

(81.0/ 

-77.2) 

85.46 

(85.2/ 

-7B.6) 

151.79 

(90.0) 

224.56 

(90.0) 

289.20 

(90.0) 

8,0 

89.04 

(90.0) 

176.19 

(90.0) 

340.18 

(90.0) 

509.10 

(90.0) 

1154.80 

(90.0) 

•  A  single  angle  applies  to  all  layers. 
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tabu  2:  Hax.»u-  fundamental  e.,enfre<t.,ncy  and  optimum  fiber  orient,, 
tions  in  degrees  of  a  four-layered  angle-ply  plate  made  of 
•Mterlal  M,  with  .  0.25,  k.|.2,3.4. 


H 

5 

10 

a 

b 

0.0 

8.26 

(0.0)* 

11.67 

(0.0) 

0.5 

8.40 

(0.0) 

11.83 

(0.0) 

I.O 

10.12 

(<5.0/ 

•45.0) 

M.83 

(A5.0/ 

-45.0) 

2.0 

19.59 

(87.8/ 

•«6.0) 

33.61 

(90.0) 

4.0 

40.29 

(90.0) 

76.86 

(90.0) 

8.0 

81.02 

(90.0) 

160.08 

(90.0) 

A  single  angle  applies  to  all 


20  <0 


13.46 

(0.0) 

14.05 

(0.0) 

14.26 

(0.0) 

13.64 

(0.0) 

14.24 

(0.0) 

14.46 

(0.0) 

17.55 

(45.0/ 

-45.0) 

18.50 

(45.0/ 

-45.0) 

18.86 

(45.0/ 

-45.0) 

47.32 

(90.0) 

54.56 

(90.0) 

57.83 

(90.0) 

132.72 

(90.0) 

187.27 

(90.0) 

223.94 

(90.0) 

.306.14 

(90.0) 

529.38 

(90.0) 

913.57 

(90.0) 

i 

»: 


r 


/ 
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TABLE  3:  HaxIauR  fundaiaefital  elgenfre<niencleSt  optlwuai  flb«roHent»tfons 
and  2(3)  ■  •iCl)  coordinates  of  a  four  layered  angle-ply  plate 
made  of  Mterlal  I 


J  5  10  20  40 


8 

H 


0.0 

9.88 

14.19 

** 

** 

** 

(Z5.9/ 

(16.0/ 

-17.7, 

-10.2, 

0.38) 

0.39) 

1.0 

13.00 

19  59 

23.73 

25.25 

25.82 

(45.0/ 

(45.0/ 

(45.0/ 

(45.0/ 

(45.0/ 

-45. C, 

-45.0, 

-45.0, 

-45.0, 

-45.0, 

0.35) 

0.35) 

0.35) 

0.35) 

0.35) 

2.0 

22.80 

39.52 

56.77 

*« 

** 

(59.4/ 

(64.1/ 

(74.0/ 

-73.3, 

‘72.3, 

-79.8, 

0.39) 

0.38) 

0.39) 

< 

4.0 

44.72 

85.77 

«* 

** 

(64.3/ 

(67.6/ 

-85.2, 

-84.7, 

0.45) 

0.44) 

The  values  are  the 

same  as 

those  given  in  Tabic  T  and  the  2(3)  ^  -  2(1) 

coordinate  has  no  effect  on  the  value  of  (u  . 


\ 


i 

f 

t 
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table  5:  HaxJmu*  separations  between  first  and  second 

second  and  third  (o.^  -  o)^)  order  elgenfrequencles  of^a  four- 
layered  angle-ply  of  mterfal  I.  with  a/h  .  lo.O. 


Objective 

“ll 

“2 

"3 

“2  -“11 

“3  •  “2 

r  •  1 

“?-“u 

18.53 

35.06 

(1.2)* 

35.06 

(2.1) 

16.53 

0.0 

(  (»i  - 

45.0, 

•  -  45.0) 

“3  -  "2 

M.M 

22.86 

(1.2) 

37.35 

(1.3) 

(2.1) 

8.12 

14.46 

(0,  - 

0.0,  » 

-  55.8) 

r  «  2 

“  *^11 

I9.ai 

40.90 

(Z.1) 

43.94 

(1.2) 

21.70 

3.04 

{6j  =  0 

'2  .0.0) 

"3  -  “2 

29.97 

45.98 

(2.1) 

67.04 

(3.1) 

(1.2 

16.02 

21.06 

(9,  •  aj.?,  0^  ,  -  55. a) 


*  The  nunbers  in  parentheses  denote  (m.n). 
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LIST  OF  FlGlIRtS 

Flg.l  Efficiency  curves  plotted  against  the  aspect  ratio  for  laminates 
made  of  material  I,  with  K=4. 

Fi9.2  Efficiency  curves  plotted  against  the  aspect  ratio  for  laminates 
made  Of  material  II,  with  K=4. 

ng.3  Percent  decrease  in  maximum  due  to  shear  deformation  plotted 
against  the  aspect  ratio,  with  K«4. 

Fig.4  Difference  in  the  efficiencies  of  four-layered  laminates  of 
material  I  optimized  with  respect  to  and  (o  ,8,). 

Fig.5  Efficiency  curves  plotted  against  the  side-to-thickness  ratio 
for  laminates  with  two,  four  and  six  layers. 

Fig.6  Efficiency  curves  plotted  against  E^/E^  ratio  for  laminates  of 
material  I,  With  K=4, 
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100(1  -Wgp/O^p) 


